Core-Level Photoelectron Spectroscopy Probing Local Strain at Silicon Surfaces and 
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Using a first-principles approach, we investigate the origin of the fine structure in Si 2p pho- 
toelectron spectra at the Si(100)2xl surface and at the Si(100)-SiO2 interface. Calculated and 
measured shifts show very good agreement for both systems. By using maximally localized Wannicr 
functions, we provide an interpretation in which the effects due to the electronegativity of second 
nearest neighbor atoms and due to the local strain field are distinguished. Hence, in combination 
with accurate modeling, photoelectron spectroscopy can provide a direct measure of the strain field 
at the atomic scale. 

PACS numbers: 79.60.-i,73.20.-r,68.35.-p 
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The increasing availability of synchrotron radiation fa- 
cilities is bearing X-ray photoemission spectra of un- 
precedented resolution characterizing surfaces and inter- 
faces The achieved sensitivity is sufficient to dis- 
tinguish inequivalent subsurface atoms with identically 
composed first-neighbor shells. Hence, the interpretation 
of such core-level spectra can no longer be achieved with 
simple electronegativity arguments, but requires the con- 
sideration of the interplay between local strain fields and 
electronegativity effects of second-nearest neighbors. 

These difficulties are strikingly illustrated for the 
Si(100)2xl surface and for the technologically relevant 
Si(100)-SiC>2 interface, which have been the objects of 
numerous highly resolved X-ray photoemission investi- 
gations. While the shifts pertaining to the first-layer 
dimer atoms of the Si(100)2xl surface have been iden- 
tified, the other lines appearing in highly resolved Si 2p 
spectra still lack a consensual assignment @, H, 0] • For 
the Si(100)-SiC>2 interface, highly resolved spectra show 
fine structure in the nonoxidized Si line, with extra com- 
ponents at lower (Si a ) and higher binding energy (SF) 
with respect to the Si bulk line 0| . 

We investigate the origin of the fine structure in Si 
2p photoemission spectra at silicon surfaces and inter- 
faces using density functional theory calculations based 
on pseudopotentials (PPs) and plane waves 0] ■ To inter- 
pret the Si 2p photoelectron spectra at silicon surfaces 
and interfaces, it is necessary to evaluate core-level shifts 
with respect to the Si bulk line. Since these shifts mainly 
result from the relaxation of valence electron states, their 
accurate determination is possible within a PP scheme, 
which does not treat core electrons explicitly [H, Q . We 
calculated Si 2p shifts including the effect of core-hole re- 
laxation by taking total energy differences between two 



TABLE I: Comparison between calculated and measured Si 
2p shifts at the Si(001)-c(4x2) surface. The shifts (in eV) are 
given with respect the Si bulk line. The meaning of labels is 
explained in Fig. [1] 
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FIG. 1: The adopted labeling of Si atoms at the Si(001)- 
c(4x2) surface. The 4' atoms can be further separated into 
inequivalent up (4'u) and down (4'd) atoms. 



separate self-consistent calculations. First, the ground- 
state energy is determined; then the PP of a given Si 
atom is replaced by another PP which simulates the pres- 
ence of a screened 2p hole in its core @|. 

For the Si(100)-c(4x2) surface we find binding energies 
of —0.49 eV and 0.02 eV for dimer-up and dimer-down 
atoms, in good agreement with experimental results 0,0, 
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FIG. 2: Calculated Si 2p shifts at the Si(100)-SiO 2 interface 
along an orthogonal direction to the interface. The Si +n cor- 
respond to oxidized Si atoms with n O neighbors (n=l - T, 
2 - A, 3 - 4 - ^) and are color-coded in the atomistic 
model. Nonoxidized Si atoms (Si ) with (o) and without (•) 
second- neighbor O atoms are distinguished. Inset: simulated 
spectrum for Si (Gaussian broadened: a— 0.08 eV) and its 
decomposition in three components. 



TABLE II: Comparison between calculated and measured Si 
2p shifts at the Si(100)-SiO 2 interface. The Si bulk line is 
taken as reference. 
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4] . The second-layer atoms are found to yield very small 
shifts (—0.01 eV) while third- and fourth-layer atoms give 
shifts to both lower (3, 4) and higher binding energies (3', 
4') with respect to the Si bulk line (Table U). Shifts of 
deeper layers are negligible. 

For the Si(100)-SiC>2 interface, we adopt a model struc- 
ture which incorporates a realistic transition region [Io| . 
Calculated binding energies are given in Fig. [2] along the 
direction orthogonal to the interface plane. The binding 
energies of oxidized Si atoms increase almost linearly with 
oxidation state Q , and show quantitative agreement with 
experimental values (Table HI)) . The focus of the present 
investigation is on nonoxidized Si atoms, which give shifts 
with a significant spread near the interface. The simu- 
lated spectrum associated to these Si atoms could be de- 
composed into three components (Fig. [21 inset), yielding 



shifts for the Si Q and SI 3 lines of -0.21 eV and 0.32 eV, 
which also agree with experimental data (Table [TT]) . 

Binding energy shifts of Si atoms with identically com- 
posed first-neighbor shells are caused by electron density 
displacements due to either the local strain fields or elec- 
tronegativity effects. We carried out an analysis in terms 
of maximally localized Wannier functions (MLWFs) 111! 
in order to distinguish between these two effects [12| . 
For the Si(100)-c(4x2) surface, we find that, apart from 
the shifts associated with the dimer atoms, the other 
shifts mainly result from the local strain induced by the 
surface reconstruction. For the interface, the Si atoms 
with second-neighbor O atoms contribute to the Si a line, 
while the Si 13 line originates from bond elongations of Si 
atoms without second-neighbor O atoms. From the ex- 
perimental shift of about 0.3 eV, we infer the occurrence 
of Si atoms with an average bond length elongation of 
- 0.05 A. 

In conclusion, we revealed the physical mechanisms un- 
derlying the fine structure in Si 2p photoelectron spectra 
at silicon surfaces and interfaces. Maximally localized 
Wannier functions offer a powerful tool to identify core- 
level shifts originating from the electronegativity of far- 
ther neighbors. A key result of our work is then that 
the remaining lines provide an atomic-scale probe of the 
strain in the structure. This confers to photoelectron 
spectroscopy a new functionality in addition to the de- 
tection of chemical composition. 
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